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Abstract. 

This paper studies two high-pollution episodes over Greater Tokyo: 9 and 10 Decem- 
ber 1999, and 31 July and 1 August 2001. Results obtained with the chemistry-transport 
model (CTM) PolairSD are compared to measurements of inorganic PM2.5. To under- 
stand to which extent the aerosol processes modeled in PolairSD impact simulated in- 
organic PM2.5, PolairSD is run with different options in the aerosol module, e.g. with/without 
heterogeneous reactions. To quantify the impact of processes outside the aerosol mod- 
ule, simulations are also done with another CTM (CMAQ). 

In the winter episode, sulfate is mostly impacted by condensation, coagulation, long- 
range transport, and deposition to a lesser extent. In the summer episode, the effect of 
long-range transport largely dominates. The impact of condensation/evaporation is dom- 
inant for ammonium, nitrate and chloride in both episodes. However, the impact of the 
thermodynamic equilibrium assumption is limited. The impact of heterogeneous reac- 
tions is large for nitrate and ammonium, and taking heterogeneous reactions into account 
appears to be crucial in predicting the peaks of nitrate and ammonium. The impact of 
deposition is the same for all inorganic PM2.5. It is small compared to the impact of other 
processes although it is not negligible. The impact of nucleation is negligible in the sum- 
mer episode, and small in the winter episode. The impact of coagulation is larger in the 
winter episode than in the summer episode, because the number of small particles is higher 
in the winter episode as a consequence of nucleation. 



1. Introduction 

With the impact of air pollution on health and vege- 
tation being a great concern, chemical transport mod- 
els (CTMs) are often used at a regional scale to pre- 
dict air quality, i.e. to compute the distribution of at- 
mospheric gases, aqueous phase species and particulate 
matter. High pollution episodes may occur depending 
on emissions and sources of pollutants, as well as on the 
dynamical characteristics of the meteorology. 

The meteorology over Greater Tokyo is strongly in- 
fluenced by orography and sea/land breeze regimes, as 
shown by Fujibe [1985] and Ohara et al. [1990]. High 
pollution episodes are often observed in early summer 
and early winter. Early winter episodes are often a con- 
sequence of the meso-synoptic scale meteorology {Mizuno 
and Kondo [1992], Uno et al. [1996], Ohara et al. [2003]): 
sea/land breeze and the blocking effect of orography 
in central Japan. These episodes are characterized 
by high NO2 and chloride concentrations ( Uno et al. 
[1996], Kaneyasu et al. [1989]). Measurements made by 
Kaneyasu et al. [1989] suggest that the precursor HNO3 
of aerosol nitrate is formed by photochemical reactions. 
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while the precursor of aerosol chloride may be locally 
emitted. For summer episodes, the role of sea breeze 
penetration for oxidants is detailed by Wakamatsu et al. 
[1999]. Hay ami [2003] founds that the PM mass com- 
pletely changes in proportion to RH in an episode that 
happened within 24 hours. Rapidly increasing RH may 
enhance condensation onto aerosols. 

Two episodes are studied in this paper: one early win- 
ter episode (9 and 10 December 1999), and one summer 
episode (31 July and 1 August 2001). These episodes 
are simulated using the CTM PolairSD (Boutahar et al. 
[2004]). This paper aims at studying the impact of dif- 
ferent aerosol processes that may influence the particle 
concentrations observed during these episodes. Further- 
more, the impact of the aerosol processes is compared 
to the impact of numerical choices made in the aerosol 
model. For example, the impact of the size distribution 
is studied by replacing the modal aerosol model MAM 
{Sartelet et al. [2006]) used in PolairSD by a sectional 
model (SIREAM: Debry et al. [2007]). Aerosol concen- 
trations are not only influenced by the aerosol model, but 
also by the parameterizations and the numerical schemes 
used for advection, diffusion, chemical mechanism {Mal- 
let and Sportisse [2006]). To have an estimate of how 
these processes may influence the concentrations, the im- 
pact of the different aerosol processes is compared to the 
impact of using a different CTM (CMAQ, Binkowski and 
Roselle [2003], Eder and Yu [2006], Yu et al. [2006b], 
Eder et al. [2006]). 

2. The Models 

The CTM Polair3D (Boutahar et al. [2004]) is used 
with the chemical mechanism RACM [Stockwell et al. 
[1997]). Photolysis rates arc computed off-line, as done 
in the photolysis rate preprocessor of CMAQ {Roselle 
et al. [1999]). Vertical diffusion is computed using the 
Troen and Mahrt's parameterization {TYoen and Mahrt 
[1986]) within the boundary layer, and using the Louis' 
parameterization {Louis [1979]) above it. PolairSD may 
be used with two aerosol models : MAM {Sartelet et al. 
[2006]) and SIREAM {Debry et al. [2007]). The difference 
between the models lies in the size distribution: in MAM 
the size distribution is modeled with four log-normal 
modes (modes i, j, k, c), and in SIREAM it is modeled 
with sections. Four sections are used in the simulations of 
this paper. In MAM and SIREAM, the modes/sections 
arc bounded as follow: mode i: < O.Ol/nm, mode j: 
[0.01; 0.1/im], mode k: [0.1; 2.5/im], mode c: > 2.5/im. 
A complete technical description of MAM and SIREAM 
may be found in Sportisse et al. [2006]. MAM is used for 
the simulations of this paper, except when specified. It 
is now briefly described. 

2.1. Composition 

Particles can be made of inert species (dust and ele- 
mental carbon), liquid water, inorganic species (sodium, 
chloride, ammonium, nitrate and sulfate) and organic 
species. 

2.2. Aerosol Processes 

Coagulation, condensation/evaporation and nuclc- 
ation are modeled as described in Sartelet et al. [2006]. 
2.2.1. Condensation/evaporation 

Condensation/evaporation is computed using the ther- 
modynamic module ISORROPIA {Nenes et al. [1998]). 
By default, thermodynamic equilibrium is assumed be- 
tween the g£is and the aerosol phases. However, Po- 
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lairSD may also be used without the assumption of ther- 
modynamic equilibrium for large modes/bins, for which 
condensation/evaporation is then computed dynamically. 
Although the thermodynamic equilibrium assumption 
may not be valid especially for large aerosols, it is used in 
the baseline simulations here because it is computation- 
ally efficient. The effect of this assumption is studied in 
a separate sensitivity simulation 

2.2.2. Nucleation 

Nuclcation is modeled using the ternary parameteri- 
zation (water, ammonium and sulfate) of Napari {Napari 
et al. [2002]). The ternary nucleation rate is several order 
of magnitudes larger than commonly used binary nucle- 
ation rates (water and sulfate). 

2.2.3. Heterogeneous reactions 

In the first runs presented here, heterogeneous reac- 
tions are not taken into account in PolairSD. If heteroge- 
neous reactions are taken into account in PolairSD, they 
are modeled according to Jacob [2000]: HO2 — > 0.5 H2O2, 
NO2 -» 0.5 HONO + HNO3, NO3 -> HNO3, N2O5 ^ 2 
HNO3. The kinetic rates of these first order reactions are 



ki = j|- ) Sa, where a is the particle radius, c? 



is the gas-phase thermal velocity in the air, Sa is the avail- 
able surface of condensed matter per air volume, and 7 
is the reaction probability that a molecule impacting the 
aerosol surface undergoes reaction. 7 strongly depends 
on the chemical and size distribution of particles. The 
values used in this paper are 7HO2 = 0.2, ^no^ = 10"'', 
7JVO3 = 10~^ and 7^205 = 0.03. 

2.2.4. Dry deposition 

Dry deposition is parameterized with a resistance ap- 
proach following Zhang et al. [2001]. The processes 
modeled include gravitational settling and the deposition 
processes of Brownian diffusion, impaction, interception, 
particle rebound. The aerodynamic resistance and the 
friction velocity are computed as in CMAQ [Bmkowski 
and Roselle [2003]). Dry deposition depends on the di- 
ameter of the particles. In Polair3D-MAM, for each log- 
normal mode and each moment, the deposition velocity is 
integrated over diameters by fourth order Gauss-Hermite 
quadratures. In PoIair3D-SIREAM, the mean diameter 
of each section is used to compute the deposition veloci- 
ties. 

2.2.5. Mode Merging versus Mode splitting 

A mode merging scheme or a mode splitting scheme is 
required in a modal model to prevent modes from over- 
lapping, i.e. to force modes to be of distinct size ranges 
throughout the simulations. Different mode merging 
schemes may be used, often based on that of Binkowski 
and Roselle [2003], where the threshold diameter between 
the two modes to be merged is chosen as the diameter 
where the number distributions of the two modes over- 
lap. In Polair3D, mode merging is applied between modes 
i and j (and between modes j and fe), when the diam- 
eter of the volume distribution of mode i (and mode k) 
exceeds a fixed diameter of O.Ol/xm for mode i (and of 
O.l/xm for mode j). 

In the tests performed later in this paper, the mode 
merging scheme is replaced by a mode splitting scheme 
{Sartelet et al. [2006]), which is designed to reproduce 
the evaluation of a mode that splits into two modes un- 
der the combined effect of nucleation, condensation and 
coagulation. 

For the simulations done with CMAQ {Binkowski and 
Roselle [2003], Mebust et al. [2003], CMAQ [1999]), the 
version 4.3 is used, with the Carbon Bond IV chemi- 
cal mechanism {Gery et al. [1989]). It is modified to 
include sodium and chloride in the Aitken and in the ac- 
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cumulation modes. In CMAQ, the aerosol module is a 
modal one with 3 modes. As in PolairSD, the following 
processes are taken into account: coagulation, conden- 
sation/evaporation, nucleation, heterogeneous reactions, 
dry deposition. For condensation/evaporation, thermo- 
dynamic equilibrium is assumed between the gas and the 
aerosol phases. Nucleation is modeled with a binary nu- 
cleation rate. The heterogeneous reaction N2O5 —^ 2 
HNO3 is taken into account. 

3. Domain and Input Data 

Simulations are performed over a 210krn x 240km area, 
centered around Tokyo, with a 5km x 5km resolution 
(Figure 1 shows the domain of simulation, which is dis- 
cretized with 42 x 48 points). The horizontal domain 
is the same for the simulations done with PolairSD and 
CMAQ. For simulations with PolairSD, 12 vertical lay- 
ers are considered (Om, 29m, 58m, 103m, 147m, 296m, 
447m, 677m, 954m, 1282m, 1705m, 2193m, 2761m). A 
no-flux boundary condition (free atmosphere) is used at 
the top boundary for diffusion. In CMAQ, 16 layers of 
varying thickness extend to about 16km height. The ver- 
tical coordinate is not altitude but sigma-lovcls. The av- 
eraged altitudes of the first 10 layers of CMAQ corre- 
spond to the altitudes used in PolairSD. Meteorological 
data axe provided by the Japanese Meteorological Agency 
with a 20km x 20km resolution every six hours. Finer 
hourly meteorological data, with a 5km x 5km resolu- 
tion are obtained by running the meteorological model 
MM5, the Fifth- Generation Pennsylvania State Univer- 
sity/National Center for Atmospheric Research (NCAR) 
Mesoscale model {Grell et al. [1994]). Initial and bound- 
ary conditions (with inputs varying every three hours) 
are obtained by running the CTM CMAQ over East 
Asia with a 45km x 45km resolution. Emission invento- 
ries are provided by a collaboration with the Japanese 
National Institute for Environmental Studies {Hayami 
and Kobayashi [2004]). Emission sources include mobile 
sources (road, air, vessels), stationary sources (domes- 
tic, industries), waste water treatment, biogenic/natural 
sources (agriculture, soil, volcanoes). Table 1 summa- 
rizes the total amount of emitted particulate matter and 
precursors over the domain. The emission inventory does 
not contain information about either the size distribution 
or the chemical speciation. The size distribution and 
the chemical speciation of PMio and PM2.5 are speci- 
fied as in CMAQ {Binkowski and Roselle [2003]). All 
PMu) — PM2.5 arc assigned to the coarse mode and par- 
ticles arc assumed to be made of 90% dust and 10% ele- 
mentary carbon. Most of PM2.5 (99.9%) are assigned to 
the accumulation mode, and 0.1% to the Aitken mode. 
For PM2.6, particles are assumed to be made of S0% dust 
and 70% elementary carbon. The parameters of the three 
modes used for emission axe shown in Table 2. Whereas 
the modal model in CMAQ has three modes (aitken, ac- 
cumulation and coarse) , the modal model in MAM has an 
additional mode designed for nanometer nucleated parti- 
cles. This mode is not used for emission. 

For gaseous species, dry deposition velocities are com- 
puted off-line following Wesely [1989]. 

Simulations start one day before the episode to allow 
for spin-up. Tests using two days for spin-up showed that 
one day spin-up is sufficient for the gaseous and PM2.5 
concentrations considered here. 
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4. Description of The Episodes 

4.1. 9*'' and 10*'* December 1999 

During this episode, chemical concentrations are ob- 
served to be high. On the 9*'', this is mostly due, to the 
presence of a meso-front. Strong cold winds coming from 
West and North West stay weak at low altitudes because 
of the presence of orographj^ in the West and in the North 
West. Warm winds from South do not penetrate much 
in land and a meso-front is observed (black dashed line 
in Figure 2). As shown in Figure 2, the meso-front is not 
so well reproduced by MM5. On the 9*'', pollutants ac- 
cumulate in the Northern part of the front. This is seen 
in Figure 3, which shows high nitrate concentrations in 
the Northern part of the front at 6pm. Sulfate concen- 
trations are also high where ammonium and nitrate con- 
centrations are high, but sulfate concentrations are even 
higher over sea. These high concentrations are brought 
from South West, but they do not propagate much over 
land because of the meso front. On the 10*'', high con- 
centrations are a consequence of weak winds and high 
pollution. 

4.2. 31** July and 1** August 2001 

The episode of 31 July and 1 August 2001 is charac- 
teristic of sea and land breeze circulation. Wind patterns 
are also influenced by local winds. 

On 31 July, winds are weak over the whole domain. 
Winds from west and south west strengthen during the 
day, and penetrate on land as the temperature on land 
gets warmer than the temperature on sea. Figure 4 dis- 
plays the wind vector at 2pm on 31 July. High sulfate 
concentrations are brought from the south west. At 2pm, 
sulfate concentrations are high over almost the whole do- 
main, as can be seen in Figure 5. 

On 31 July at night, although south westerly winds 
are strong over sea, winds stay low on land. In the early 
morning of 1 August, south westerly winds weaken. They 
do not penetrate much on land, and winds on land stay 
low. 

5. Comparison to Measurements 
5.1. Measurements 

Aerosol distributions were collected using impactor, 
denuder, and filter pack [Hayami and Fujita [2004]), 
every three hours, except for the summer 2001 where 
at night they were collected every six hours. Hourly 
gaseous distributions at different stations distributed over 
Greater Tokyo are provided by Japanese local institute 
governments. 

The stations at which comparisons are made are shown 
in Figure 1. 

On 9 and 10 December 1999, for gas, data are available 
for O3, NOx and SO2 at 3 sites: Yokosuka, Fukaya and 
Kudau. For PM2.5, data arc available for sulfate, ammo- 
nium, nitrate and chloride at 4 sites: Yokosuka, Omiya, 
Fukaya and Kudan (Tokyo). 

On 31 July and 1 August 2001, data are available at 
two sites: Kumagaya and Komae, for O3, NOa,, SO2 and 
PM2.5. 

Sulfate concentrations are overall higher on 31 July 
and 1 August 2001 than on 9 and 10 December 1999, 
with an average over the stations of 13.3/xg m~^ against 
2.6/ig m~^. Furthermore, sulfate represents as much 
as 50% of the inorganic PM2.5 in the summer episode, 
against 21% in the winter one. Ammonium and nitrate 
concentrations are also higher on 31 July and 1 August 
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2001 than on 9 and 10 December 1999, with an aver- 
age of 7.7/jg against 2.7/Lig for ammonium, and 
5.7/ig m~^ against 4.5/ig for nitrate. 

5.2. Statistics 

Comparison of the results obtained with PolairSD to 
measurements is done for inorganic fine aerosols (PM2.5). 
Model results for this comparison take into account data 
only from the grid boxes for which observations are 
available. In MAM/SIREAM, PM2.5 are computed by 
summing the three smallest modes/sections (modes i, j 
and k). Similarly, in CMAQ, as the mode i is not mod- 
eled, the two smallest modes, i.e. j and k, are summed 
up to compute PM2.5 

The relative bias and error between modeled and ob- 
served concentrations arc quantified using unbiased sym- 
metric metrics ( Yu et al. [2006] ) : the normalized mean 
bias factor Batmsf and the normalized mean absolute 
error factor 'Enmaef (Appendix A). The correlation co- 
efficient (%) is also used as a statistical indicator. The 
smaller the bias and the error are, and the larger the 
correlation is, the closer the model fits the observation. 
Bias indicates whether the model tends to under or over- 
predict the observation, and error indicates how large the 
deviation is. 

Yu et al. [2006] suggests the criteria of model perfor- 
mance for sulfate to be taken as \Bnmbf\ < 25% and 
Enmaef < 35%. 

5.3. 9*^ and 10*'' December 1999 

The statistical indicators are shown in Table 3. The 
model satisfies the criteria suggested by Yu et al. [2006] 
for sulfate. However, these criteria are not satisfied for 
ammonium, nitrate and chloride, which are more diffi- 
cult to model because of their volatility. This difficulty 
to model nitrate and ammonium is stressed for example 
by Zhang et al. [2006] where nitrate and ammonium are 
undcrprcdictcd by factors 9.6 and 2.1 respectively in the 
Southeastern US for the period of 1-10 July 1999. 

The correlation coefficients are good for all species, 
ranging from 66% for sulfate to 36% for chloride. This 
suggests that the overall diurnal patterns are well mod- 
eled. As shown in Figure 8, the majority of hourly sim- 
ulations falls within a factor 2 of the observations for 
sulfate and for high values of nitrate and ammonium. 
However, values lower than about 3/ig m~^ for ammo- 
nium and 5^g m^"* for nitrate arc often scattered outside 
the factor of two reference line. 

The presence of the meso-front in the afternoon of 9 
December may be seen from the concentrations of PM2.5 
in Figures 6 and 7. Nitrate and ammonium concentra- 
tions axe high in the afternoon before 6pm at Omiya and 
Fukaya, which axe located in the northern part of the 
front where the pollutants accumulate. However, they 
are low at Yokosuka, which is located in the south part 
of the front. Kudan is located close to the front at 6pm. 
Concentrations arc low at Kudan before 6pm and the 
concentrations of pollutants increase at Kudan from 6pm. 
On 10 December, pollutants are observed to be high at 
the four stations, as a consequence of weak winds. 

PolairSD tends to over-estimate sulfate as shown by 
the Bnmbf, which is as high as 0.26. As shown in Fig- 
ures 6 and 7, sulfate is overestimated in the evening of 
the 10*''. For both ammonium and nitrate, although the 
position of the peaks in time is well predicted by Po- 
lairSD, they arc under-estimated at Fukaya and Omiya 
on the 9*''. To find the reason of the discrepancies be- 
tween simulations and measurements on the 9*'' for the 
peaks of nitrate and ammonium at Fukaya and Omiya, 
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the sensitivity of the amphtude of these peaks to op- 
tions in the aerosol module is studied in the next section. 
Ammonium is overall over-estimated, while nitrate and 
chloride are under-estimated. However, for both nitrate 
and ammonium, the model has difHculties to predict the 
low concentrations in the night between the 9"" and the 
10*'' at Omiya and Fukaya, as confirmed by the scatter 
plots. The reason for this discrepancy is also investigated 
in the next section. 

5.4. 31** July and 1** August 2001 

Comparisons of the results obtained with PolairSD are 
shown in Figure 9 for PM2.5. 

The statistical indicators arc shown in Table 4 for 31 
July and 1 August 2001. As for 9 and 10 December 1999, 
the model satisfies the criteria suggested by Yu et al. 
[2006] for sulfate. However, these criteria are again not 
satisfied for ammonium and nitrate. Because chloride 
concentrations are very low in this summer run, they are 
not shown. 

As for the winter episode, the majority of hourly sim- 
ulations falls within a factor 2 of the observations for sul- 
fate (Figure 8). The majority of hourly simulation falls 
within a factor 2 of the observation at Komac, but scat- 
ter outside the reference lines is observed at Kumagaya. 
Considerable scatter is observed for nitrate at both Ko- 
mae and Kumagaya with many simulated values falling 
outside the factor of two reference lines. 

The correlation coefficient is low for sulfate (-2%). For 
example, the model predicts a decrease in sulfate concen- 
trations in the morning of 31 July at both Komae and 
Kumagaya (Figure 9), which is not confirmed by obser- 
vations. Because this episode is a clear-sky one, cloud 
chemistry does not infiuence sulfate concentrations. Sul- 
fate may be produced by the condensation of H2SO4, or 
it may be transported to the domain of study through 
boundary conditions. H2SO4 is either directly emitted 
or produced by the reaction of 502 with OH. How- 
ever, although the correlation coefficient is low for sulfate, 
good correlation is observed for SO2 (Table 4). This sug- 
gests that uncertainty in sulfate concentrations is linked 
to uncertainty in the H2SO4, emissions and in the sulfate 
boundary conditions. 

The sulfate concentrations, which are quite high all 
through the episode, are relatively well modeled. How- 
ever, the measured concentrations do not vary much, 
whereas the model exhibits stronger time variations. For 
ammonium, although results are close to measurements 
at Komac, ammonium is under-estimated at Kumagaya. 
According to measurements, ammonium concentrations 
are twice as high in Kumagaya as in Komae. These high 
concentrations in Kumagaya are not reproduced by the 
models. Nitrate is severely under-estimated at both Ku- 
magaya and Komae. The reasons for the discrepancies 
between measurements and modeled concentrations are 
investigated in the next section. 

6. Impact of Aerosol Processes 

The impact of different processes on aerosol concen- 
trations is investigated for both episodes. The following 
processes are considered: rmcloation, coagulation, con- 
densation/evaporation, dry deposition, heterogeneous re- 
actions. The impact of numerical schemes in the aerosol 
module is also studied: size-distribution and mode merg- 
ing/splitting. 

The impact of each of these physical and numerical 
processes is quantified by comparing a simulation Sp 
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where only a process p is ignored (or its parameteriza- 
tion is changed) to the reference simulation where 
all the processes are taken into account. The simula- 
tions Sp are successively: a simulation without conden- 
sation/evaporation, a hybrid simulation (thermodynamic 
equilibrium is only assumed for the smallest mode instead 
of being assumed for the four modes) , a simulation with 
heterogeneous reactions (heterogeneous reactions are not 
taken into account in the reference simulation Sr), a sim- 
ulation with only the N2O5 heterogeneous reaction, a 
sinnilation without the NO2 heterogeneous reaction, a 
simulation without deposition, a simulation without nu- 
cleation, a simulation without coagulation, a simulation 
done using SIREAM instead of MAM (the log-normal 
size distribution is replaced by a sectional distribution), 
a simulation done using mode splitting instead of mode 
merging, and a simulation without neither mode splitting 
nor mode merging. The quantification of the impact of 
each process p is done by computing the Y^mmbf and the 
'EiNMAEF between Sp and Sr at the sites where measure- 
ments are made. The larger in absolute value the Bjvmbf 
and the 'Enmaef arc, the stronger the impact of the pro- 
cess p is. The statistics Bjvmbf and ^nmaef between 
CMAQ and Sr axe also computed in order to compare the 
impact of using a different chemistry transport model to 
the impact of each of the processes p. CMAQ takes into 
account condensation/evaporation with the same ther- 
modynamic model as PolairSD (Isorropia), the heteroge- 
neous reaction of N2O5, deposition, nuclcation and mode 
merging. The comparisons of simulations done with Po- 
lairSD allow us to quantify the uncertainty related to the 
aerosol module. By using CMAQ, the simulated PM2.5 
concentrations may vary not only because of uncertainty 
in the aerosol module, but aJso because of uncertainties 
in chemistry, transport and diffusion. For example, even 
though PolairSD and CMAQ use the same meteorological 
fields, uncertainties in transport may be linked to the nu- 
merical schemes used, uncertainty in diffusion is related 
to the modeling of vertical diffusion. 

Because the reference run does not always capture the 
observed behavior or because it may capture it for wrong 
reasons (e.g. inaccuracies in several processes cancel each 
other out), changing processes one at a time does not 
automatically give a reliable estimate of the relative im- 
portance of these processes in the real world, but only in 
the model. Similarily, the comparison of two CTMs docs 
only give an estimate of the uncertainty related to trans- 
port, chemistry and diffusion, as these could be described 
inaccurately in both models. 

Tables 5 and 6 compare the reference simulation Sr 
to different simulations Sp for 9 and 10 December 1999 
and SI July and 1 August 2001 respectively. For both 
episodes, the impact of heterogeneous reactions is large 
for nitrate and ammonium (the 'Enmaefs are as large 
as 1.95 and 0.79 respectively in winter). The HNO3 pro- 
duced by heterogeneous reactions condenses onto aerosols 
to form nitrate, and the available ammonia condenses to 
neutralize the nitrate. Concentrations of nitrate and am- 
monium increase by taking into account heterogeneous 
reactions, as shown by the positive 'Qnmbf- The im- 
pact of heterogeneous reactions on sulfate is small. The 
H2O2 concentration produced by heterogeneous reactions 
does not influence sulfate concentrations because aque- 
ous chemistry is not taken into account. The impact 
of condensation/evaporation is preponderant for ammo- 
nium, nitrate and chloride (for nitrate, the Enmaef is 
as large as 2.44 in the winter episode and 96.4 in the 
summer episode). This impact is even larger than the 
impact of using CMAQ, or the impact of heterogeneous 
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reactions. Condensation is largely preponderant over 
evaporation as seen by the negative Bjvmbf- Because 
of the low-volatility of sulfate, the impact of the hybrid 
scheme is small for sulfate (with a Enmaef under 0.07). 
It is larger for other inorganic species, such as ammo- 
nium and nitrate, which are influenced by whether con- 
densation/evaporation is computed dynamically or by as- 
suming thermodynamic equilibrium. The hybrid scheme 
mostly influences ammonium in the winter case with a 
EjvMAEF of 0.41, and it influences mostly nitrate in the 
summer case with a Enmaef of 0.49. Because coagula- 
tion and deposition do not differentiate the chemical com- 
position of aerosols, their impacts are the same for each 
chemical species. The impact of deposition is between 
0.08 and 0.10 for each species for both the winter and 
the summer simulations. Not taking deposition into ac- 
count, leads as expected to an increase in aerosol concen- 
trations as shown by the positive Unmbf- The impact of 
nucleation is small, with Enmaef smaller than 0.06. Al- 
though a large number of nanometer particles arc created 
by nucleation (the ternary nucleation scheme is several 
order of magnitudes larger than commonly used binary 
nucleation schemes), the mass produced is small com- 
pared to the mass of PM2.6. Furthermore, although the 
nucleated particles are made of sulfate and ammonium, 
nitrate and chloride are also influenced by nucleation, 
as they may condense onto freshly nucleated particles. 
Considering numerical processes, the impact of mode 
splitting and mode merging is negligible with Enmaef 
smaller than 0.03. Mode merging/splitting influence the 
smaller particles, which do not contribute much to the 
mass of PM2.5. The impact of using SIREAM, although 
not negligible, is not large for sulfate with Enmaef un- 
der 0.09. The impact is larger for nitrate, especially in 
the summer episode, where the Enmaef reaches 0.22. A 
small impact on sulfate concentrations may correspond 
to larger impacts on other inorganic scmi-volatilc com- 
ponents. During condensation/evaporation, ammonium, 
nitrate and chloride arc not only influenced by the size 
distribution of aerosols, but also by the internal compo- 
sition because of their volatility. 

6.1. 9*^ and 10*'' December 1999 

In the winter episode of 1999, for sulfate, there is a 
dominant impact of both condensation and coagulation 
with Enmaef as large as 0.27 and 0.21 respectively. 
These impacts are of the same order of magnitude as 
those of CMAQ (the Enmaef is 0.23). For sulfate, the 
impacts of other processes are small, except for deposi- 
tion, for which Enmaef is equal to 0.10. Sulfate is pro- 
duced by the condensation of H2SO4 and by the transport 
from outside the model domain through boundary con- 
ditions. Because sulfate decreases only by a factor 1.27 
if condensation is not taken into account, the impact of 
long-range transport is likely to be high. The boundary 
condition for sulfate averaged over time, latitude and lon- 
gitude is 3.0 /Ltg m~^, against a mean concentration of 3.5 
fig m~^ in the domain of study. For other species, the 
impact of condensation/evaporation largely dominates, 
with Enmaef as largo as 3.92 for chloride, 2.44 for ni- 
trate and 1.04 for ammonium. This indicates that chlo- 
ride, nitrate and ammonium are strongly influenced by 
local conditions. The impact of heterogeneous reactions 
is important for nitrate and ammonium with a Enmaef 
as large as 1.95 and 0.79. The dominant heterogeneous 
reactions are the N2O6 and the NO2 heterogeneous reac- 
tions. The Enmaef is as high as 0.94 and 0.38 for nitrate 
and ammonium when only the N2O5 heterogeneous reac- 
tion is taken into account. A simulation without the NO2 
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heterogeneous reaction but with all the other three gives 
similar results to the simulation with only the N2O5 het- 
erogeneous reaction. The impact of the hybrid scheme is 
important for ammonium with a 'Enmaef of 0.41, which 
is as large as the impact of using CMAQ (0.38). However, 
the impact is smaller for other chemical species, with a 
Eatmabf of 0.13 for nitrate against a ^nmaef of 1.0 
when CMAQ is used. Coagulation has an impact almost 
as large as the impact of condensation or the impact of 
CMAQ for sulfate with a Ejvmabf of 0.21. Although 
the impact of coagulation is the same for each chemical 
species, it is less important for nitrate, ammonium and 
chloride in comparison to the impact of other processes. 
Although nucleation is small, it is not negligible with a 
Eatmabf of 0.05 for sulfate and ammonium. 

The impacts of the dominant processes during the 
episode are now investigated in more details, aiming at 
understanding the potential reasons for the discrepancies 
between measurements and simulated concentrations at 
the different stations (Figures 6 and 7). As seen in section 
5, sulfate is overall ovcrpredicted. This is most of the time 
due to uncertainties in condensation and coagulation. For 
example, at Omiya and Kudan, Polair3D tends to over- 
estimate the sulfate concentrations but CMAQ does not. 
Simulations without condensation or without coagulation 
give results that are similar to CMAQ, although both 
condensation and coagulation arc taken into account in 
CMAQ. As discussed in section 5, the peaks of nitrate 
and ammonium are undorpredicted on the 9*'' at the sta- 
tions above the meso-front: Omiya and Fukaya. When 
heterogeneous reactions are taken into account, the peaks 
are not underestimated any longer. However, they tend 
to be predicted earlier, stressing the difficulties of MM5 
to accurately simulate the meso-front. Furthermore, the 
differences between Polair3D and CMAQ are often ex- 
plained by the N2O5 heterogeneous reaction. If this het- 
erogeneous reaction is taken into account, both PolairSD 
and CMAQ predict high nitrate and ammonium concen- 
trations in the night and morning of the 9*'' at all four 
stations and in the night and in the morning of the 10*'' 
at Kudan, Omiya and Yokosuka. However, the high con- 
centrations on the 10*** are not observed, suggesting the 
need to revise the rate of the N2OB heterogeneous re- 
action. For example, Evans and Jacob [2005] suggest 
a rate that vary with the aerosol composition, temper- 
ature and relative humidity. The low nitrate and am- 
monium concentrations between the 9"' and the lO"' at 
Omiya are ovcrpredicted by CMAQ because of the het- 
erogeneous reactions. However, although the low nitrate 
concentrations are quite well predicted by Polair3D when 
heterogeneous reactions are ignored, the low ammonium 
concentrations are still overpredicted. Although nitrate 
concentrations are very small all through the episodes if 
condensation is not taken into account, ammonium con- 
centrations are larger than the lowest observed concen- 
trations. The unability of the model to reproduce the 
low ammonium concentrations may therefore be linked 
to uncertainties outside the aerosol module, most likely 
in the boundary conditions. 

6.2. 31** July and 1** August 2001 

In the summer episode, for sulfate, the impacts of 
the different processes are small compared to the impact 
of using CMAQ. The 'Enmaef for the simulation with 
CMAQ is equal to 0.55 against 0.08 for other simula- 
tions such as the one without condensation, the one using 
SIREAM and the one without deposition. The impact of 
condensation on sulfate is not as strong as it is in the 
winter episode. This suggests that sulfate mostly comes 
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from long-range transport. Whereas in the winter case, 
the averaged sulfate in the domain is larger than the av- 
eraged sulfate used for boundary conditions, the opposite 
holds in the summer case. The boundary condition for 
sulfate averaged over time, latitude and longitude is 16.6 
fj,g m~*, against a mean concentration of 13.5 /ig m~^ in 
the domain of study. In fact, Mount Oyama of Miyake 
Island, which is located about 180km south of central 
Tokyo, was in eruption at that time [Shinohara et al. 
[2003]). Sulfate concentrations may be more sensitive to 
transport or diffusion processes, such as the parameter- 
ization used to model the vertical diffusion [Mallet and 
Sportisse [2006]) However, condensation influences to a 
great extent the nitrate and ammonium concentrations, 
with a EjvMABF as large as 96.4 and 9.4 respectively. The 
impact of condensation on ammonium and nitrate is large 
compared to the impact of using CMAQ. Ammonium and 
nitrate arc produced locally. Although condensation is 
largely a dominant process for the ammonium and ni- 
trate concentrations, the influence of the hybrid simula- 
tion, where the thermodynamic equilibrium assumption 
is removed, is limited with a Enmaef of only 0.49 for 
nitrate and 0.18 for ammonium. The relative importance 
of the hybrid scheme versus condensation is higher in the 
winter episode than in the summer episode. Perhaps, this 
is because thermodynamic equilibrium is reached quickly 
under high temperatures (the mean temperature is 295K 
in the summer episode and 278K in the winter episode) 
and high pollutant concentrations (the averaged PM2.5 
concentration is 20.1^ig m^'* in the summer episode and 
13.3/^g m~^ in the winter episode) ( Wexler and Seinfeld 
[1990]). Nitrate concentrations, and ammonium concen- 
trations to a lesser extent, axe also strongly influenced by 
heterogeneous reactions with a F,nmaef of 1.20 and 0.12 
respectively. As for the winter case, the dominant het- 
erogeneous reactions arc the N2O5 and the NO2 hetero- 
geneous reactions. The Enmaef for the N2O5 heteroge- 
neous reactions is as large as 0.95 for nitrate and 0.09 for 
ammonium. The impact of using SIREAM is not small 
although not preponderant: it is as large as the impact 
of deposition and condensation for sulfate, and it is as 
large as the impact of the hybrid scheme for ammonium 
and nitrate. Mode splitting, mode merging, nucleation 
and coagulation are negligible during this episode. Nucle- 
ation and coagulation are negligible in the summer case 
but not in the winter case. The mean temperature and 
relative humidity in the winter case are 278K and 49%, 
while they are 295K and 76% in the summer case. Ac- 
cording to Korhonen et al. [2003], under these conditions 
of temperature and relative humidity, the nucleation rate 
is about 6 times higher in the winter case than in the sum- 
mer case. The impact of coagulation is larger than the 
impact of nucleation, because coagulation does not only 
influence particles when they arc freshly nucleated, but 
coagulation also influences those particles as they grow 
by condensation. 

The impacts of the dominant processes during the 
episode are now investigated in more details, aiming at 
understanding the potential reasons for the discrepan- 
cies between measurements and simulated concentrations 
at the different stations (Figure 9). The sulfate con- 
centrations are very little sensitive to the options used 
in the aerosol module in PolairSD, but large differences 
are observed between CMAQ and Polair3D. For example, 
at Kumagaya, Polair3D predicts lower sulfate concentra- 
tions than CMAQ in the night and morning of the 31"', 
whereas it predicts higher sulfate concentrations later in 
the episode. Furthermore, the models exhibit stronger di- 
urnal variations compared to observations. Because sul- 
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fate concentrations mostly come from long-range trans- 
port, and are little sensitive to options used in the aerosol 
module, the uncertainties in the sulfate concentrations 
are linked to uncertainties in the meteorology and bound- 
ary conditions. PolairSD and CMAQ use the same mete- 
orological fields, except for the vertical diffusion, which is 
computed following TYoen and Mahrt [1986] in PolairSD 
and following CMAQ [1999] (K-theory) in CMAQ. Un- 
certainty due to vertical diffusion is likely to be important 
as stressed by Mallet and Sportisse [2006]. Boundary 
conditions are provided on the vertical grid of CMAQ, 
which is used for the continental run, and they are pro- 
jected onto the vertical grid of PolairSD. Uncertainties 
due to boundary conditions may not only be linked to this 
projection, but also to the projection of moments. Be- 
cause the boundary conditions are obtained from CMAQ, 
they arc given as the moments of order zero, two and 
three for particulate matter. However, a projection is 
required by PolairSD which uses the moments of order 
zero, three and six to represent particulate matter. Am- 
monium and nitrate are strongly influenced by conden- 
sation. Their concentrations arc very low, and almost 
zero for nitrate, if condensation is not taken into account. 
As ammonia may condense onto aerosols to neutralize 
the sulfate, ammonium follows the diurnal evolution pre- 
dicted for sulfate (the correlation between computed sul- 
fate and ammonium is as high as 87%). Errors on am- 
monium concentrations arc partly due to errors on sul- 
fate concentrations. Another cause of errors may come 
from uncertainties in total ammonium. Ammonium con- 
centrations are better predicted at Komae than at Ku- 
magaya. As shown in Figure 10, the total ammonium 
concentration computed by PolairSD is higher at Komae 
than at Kumagaya, even though ammonium concentra- 
tion in the eierosol phase is higher at Kumagaya. In the 
first part of the episode, at Kumagaya, the total ammo- 
nium computed is as low as the measured ammonium in 
the aerosol phase. Even though ammonium is system- 
atically underestimated at Kumagaya, some discrepan- 
cies between modeled and observed ammonium can be 
explained by processes in the aerosol module at Komae. 
For example, whereas PolairSD overestimates ammonium 
concentrations in the afternoon of the SI"*, the ammo- 
nium concentrations are underestimated if the hybrid 
scheme is used. Nitrate is strongly undcrprcdictcd by 
PolairSD. However, the peaks of nitrate at night are bet- 
ter predicted when the N2O5 iictcrogcneous reaction is 
taken into account, although these peaks are sometimes 
overprcdictcd by CMAQ. As for total ammonium, total 
nitrate concentration is much higher at Komae than at 
Kumagaya (Figure 11). However, even though the total 
nitrate concentration is high at Komae, the high concen- 
trations of nitrate in the aerosol phase in the afternoon of 
the 31"* at Komae are not reproduced even if the hybrid 
scheme is used. Yu et al. [2005] suggest that measure- 
ment uncertainties in sulfate and total ammonium, i.e. 
ammonium plus ammonia, can account for most of the 
discrepancies between the model predictions and obser- 
vations in partitioning of aerosol nitrate. 

7. Conclusion 

In this article, two high-pollution episodes over 
Greater Tokyo are studied: one during the winter 1999 
(9 and 10 December) and one during the summer 2001 
(SI July and 1 August 2001). For each of these episodes, 
the chemistry transport model PolairSD is compared to 
measurements for inorganic components of PM2.B. 

For sulfate, error statistics are in agreement with 
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model performance criteria {Yu et al. [2006]). Inorganic 
components of PM2.5 remain overall well modeled except 
for nitrate in the summer episode. 

To understand to which extent the aerosol processes 
modeled in PolairSD influence the particle concentrations 
during the summer and the winter episodes, different sim- 
ulations are made where only one process differs from 
the default or reference simulation. The following phys- 
ical/chemical processes arc considered: nucloation, co- 
agulation, condensation/evaporation, whether condensa- 
tion is modeled dynamically or using the thermodynamic 
equilibrium assumption, dry deposition, heterogeneous 
reactions. For numerical processes, the impact of mode 
merging/mode splitting and the impact of the size distri- 
bution (modal versus sectional) axe evaluated. A compar- 
ison of the impact of each aerosol process described above 
to the impact of using the CTM CMAQ allows us to as- 
sess the importance of using different parameterizations 
and numerical schemes not only for aerosol processes but 
also for chemistry, transport and diffusion. 

This study illustrates that the impact of aerosol pro- 
cesses on aerosol concentrations differs depending on lo- 
cal conditions and aerosol chemical components. For ex- 
ample, in the summer episode, for sulfate, the impact 
of long-range transport largely dominates. In the win- 
ter episode, sulfate is mostly impacted by condensation, 
coagulation, long-range transport, and deposition to a 
lesser extent. Whereas nucleation and coagulation are 
negligible in the summer episode, they arc not in the win- 
ter episode. The impact of coagulation is larger in the 
winter episode than in the summer episode, because the 
number of small particles is higher in the winter episode 
as a consequence of nucleation. The impact of conden- 
sation/evaporation is dominant for ammonium, nitrate 
and chloride in both episodes. However, the impact of 
the thermodynamic equilibrium assumption is limited. 
The impact of heterogeneous reactions is large for nitrate 
and ammonium. The dominant heterogeneous reactions 
are the NO2 ^ 0.5 HONO + HNO3 and the N2O5 ^ 
2 HNO3 reactions. The impact of using a sectional rep- 
resentation of the size distribution is not negligible, and 
it is higher for ammonium and nitrate than for sulfate. 
The impact of mode merging/mode splitting is negligible 
in both episodes. 

The comparison of the different runs also allows us 
to understand discrepancies between observed and sim- 
ulated inorganic PM2.5 at different stations. Heteroge- 
neous reactions appear to be crucial in predicting the 
peaks of nitrate and ammonium in the winter episode. 
However, heterogeneous reactions sometimes lead to con- 
centrations that are too high, suggesting the need for a 
more detailed parameterisation of the reaction rates. 

Although the impact of mode merging and mode split- 
ting is negligible on PM2.5 concentrations, it may influ- 
ence the size distribution of aerosols. Larger differences 
between the different runs may be observed by comparing 
the size distribution of aerosols or the concentrations of 
smaller particles such as PMi. In particular, the impact 
of using a sectional rather than a modal model would be 
larger on the size distribution than on the mass of PM2.5. 
The impact of nucleation and whether condensation is 
computed dynamically or using the thermodynamic as- 
sumption would be larger on PMi than on PM2.5. 
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Appendix A: Statistical Indicators 

The following indicators are computed in order to eval- 
uate error statistics for model-to-data comparisons. Let 
(oi)^ and (cf)j be the observed and the modeled concen- 
trations, where i is over n time series and locations, and 
o = Y^^=i °i ^ ~ Yli=i averaged observed and 

modeled concentrations respectively 

• Correlation 



Er=l - C) (Oi - O) 



(Al) 



\/Er=i - c)' \/Er=i {oi-of 
1 ^ 1 ^ 

with: o = — > Oi and c = — > c, (A2) 
n n ^ 



n ^-^ n ■ 

i=l i=l 



Normalized mean bias factor {Bnmbf) 



3 - 1 if c > o (A3) 
o 

1 - 3 if c < o (A4) 
c 



• Normalized mean absolute error factor {E,nmaef) 



o 



if c > o (A5) 
if c < o (A6) 
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Figure 1. Location of stations at which the comparisons 
to data are made. 



Table 1 . Total amount (in kt) of emitted particulate matter 
and precursors over the domain (210km x 240km) for each 
episode. 





PM2.5 


PMio 




50x 


NH3 


HCl 


Winter episode 
Summer episode 


0.44 
0.40 


0.23 
0.23 


3.72 
2.82 


1.75 
1.24 


0.54 
1.81 


0.16 
0.15 
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Table 2. Parameters of the three modes used for emission. 





Aitken 


Accumulation 


Coarse 


Mean Diameter (//m) 


0.03 


0.3 


6 


Standard Deviation 


1.7 


2 


2.2 
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Figure 6. PM2.5 concentration of sulfate and nitrate at 
Fukaya, Yokosuka, Kudan, Omiya for 9 and 10 December 
1999. 
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Figure 7. PM2.5 concentration of ammonium and chlo- 
ride at Fukaya, Yokosuka, Kudan, Omiya for 9 and 10 
December 1999. 
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Table 3. Correlation (corr), Bjvmbf and Ejstmaef ob- 
tained with PolairSD for 9 and 10 December 1999. 





corr 


^NMBF 


^NMAEF 


Sulfate 
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Chloride 
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Figure 8. Scatter plots of observation (ordinate) versus 
simulation (abscissa) for sulfate, ammonium and nitrate 
for 9 and 10 December 1999 (left panels) and for 31 July 
and 1 August 2001 (right panels). 1:1, 1:2, and 2:1 refer- 
ence lines are provided. 
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Figure 9. Comparisons of the PM2.5 concentrations of 

sulfate, nitrate and ammonium at Kumagaya and Komae 
for 31 July and 1 August 2001 for different simulations. 
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Table 4. Correlation (corr), Bjvmbf and Ejvmabf ob- 
tained with PolairSD and CMAQ for 31 July and 1 August 
2001. 
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Table 5. Comparison of Polair3D-MAM to different simula- 
tions for 9 and 10 December 1999: Bjvmsf (B) and Enmaef 
(E). 
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Table 6. Comparison of Polair3D-MAM to different sim- 
ulations for 31 July and 1 August 2001: ^nmbf (B) and 



Sulfate Ammonium Nitrate 
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Figure 10. Concentrations of ammonium (modeled and 

observed) and total ammonium at Kumagaya and Komae 
for 31 July and 1 August 2001. 
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Figure 11. Concentrations of nitrate (modeled and ob- 
served) and total nitrate at Kumagaya and Komae for 31 
July and 1 August 2001. 
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